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Nomenclature 


ACE 

= Advanced Composition Explorer 

D 

= Dose (cGy) 

E 

= Kinetic energy (MeV) 

GCR 

= Galactic Cosmic Ray 

GeV 

= Giga electron Volt 

H 

= Hydrogen 

HZETRN 

= High charge (Z) and Energy TRaNsport code 

ICES 

= International Conference on Environmental Systems 

ISS 

= International Space Station 

LaRC 

= Langley Research Center 

LEO 

= Low Earth Orbit 

LIS 

= Local Interplanetary Spectrum 

MC 

= Monte Carlo 

MeV 

= Mega electron Volt 

Ni 

= Nickel 

NIST 

= National Institute of Standards and Technology 

PE 

= Polyethylene (CH 2 ) 

PET 

= Polyethylene Terephthalate (C 10 H 8 O 4 ) 

PP 

= Polypropylene 

RBE 

= Relative Biological Effectiveness 

R(50) 

= Range at 50 MeV protons (g/cm 2 ) 

S 

= Stopping power (MeV/g/cirf) 

SAE 

= Society of Automotive Engineers 

SAMPE 

= Society for the Advancement of Material and Process Engineering 

SPE 

= Solar particle event 

Z 

= Particle charge number 

9 

= Differential GCR flux [#/(MeV-cm 2 -day)] 
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Abstract 


The origin of the aluminum equivalent shield approximation in space 
radiation analysis can be traced back to its roots in the early years of the NASA 
space programs (Mercury, Gemini and Apollo) wherein the primary 
radiobiological concern was the intense sources of ionizing radiation causing 
short term effects which was thought to jeopardize the safety of the crew and 
hence the mission. Herein, it is shown that the aluminum equivalent shield 
approximation, although reasonably well suited for that time period and to the 
application for which it was developed, is of questionable usefulness to the 
radiobiological concerns of routine space operations of the 21 st century which 
will include long stays onboard the International Space Station (ISS) and 
perhaps the moon. This is especially true for a risk based protection system, as 
appears imminent for deep space exploration where the long-term effects of 
Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis 
demonstrates that sufficiently large errors in the interior particle environment of 
a spacecraft result from the use of the aluminum equivalent approximation, and 
such approximations should be avoided in future astronaut risk estimates. In this 
study, the aluminum equivalent approximation is evaluated as a means for 
estimating the particle environment within a spacecraft structure induced by the 
GCR radiation field. For comparison, the two extremes of the GCR environment, 
the 1977 solar minimum and the 2001 solar maximum, are considered. These 
environments are coupled to the Langley Research Center (LaRC) deterministic 
ionized particle transport code High charge (Z) and Energy TRaNsport 
(HZETRN), which propagates the GCR spectra for elements with charges (Z) in 
the range 1 <Z <28 (H - Ni) and secondary neutrons through selected target 
materials. The coupling of the GCR extremes to HZETRN allows for the 
examination of the induced environment within the interior of an idealized 
spacecraft as approximated by a spherical shell shield, and the effects of the 
aluminum equivalent approximation for a good polymeric shield material such as 
generic polyethylene (PE). The shield thickness is represented by a 25 g/cm ' 
spherical shell. Although one could imagine the progression to greater thickness, 
the current range will be sufficient to evaluate the qualitative usefulness of the 
aluminum equivalent approximation. Upon establishing the inaccuracies of the 
aluminum equivalent approximation through numerical simulations of the GCR 
radiation field attenuation for PE and aluminum equivalent PE spherical shells, 
we further present results for a limited set of commercially available, hydrogen 
rich, multifunctional polymeric constituents to assess the effect of the aluminum 
equivalent approximation on their radiation attenuation response as compared to 
the generic PE. 
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1. Introduction 


The space radiation concerns of the early years of the US space programs (Mercury, 
Gemini and Apollo) were for the more intense sources of ionizing radiation, requiring protection 
to avoid radiation sickness which could impact mission safety or even cause death during the 
mission (ref. 1). The Relative Biological Effectiveness (RBE) for such effects was assumed to 
be near unity, so that the energy absorbed per unit mass of tissue (i.e. dose) was the indicator of 
health prognosis, which led to considerable simplification in dosimetry analysis in which tissue 
equivalent ion chambers played a dominant role. The primary environments of concern were for 
the geomagnetically trapped radiation belts (i.e. protons and electrons) and high flux Solar 
Particle Events (SPE), in which protons are the largest constituents and contributions of 
secondary radiations to the dose from both trapped and SPE environments are limited. For 
example, a 100 MeV proton will have a few percent probability of a nuclear reaction in coming 
to rest by a cascade of atomic interactions. These considerations were carried out in the 
spacecraft shield design in which nuclear reactions were normally ignored in the design process 
(refs. 2, 3). For historical reasons, it is important to note that the exclusion of nuclear processes 
contribution to the dose to evaluate the secondary radiation field in even a simplified geometry 
was in part motivated by the very long computational time required by the available Apollo era 
statistical based Monte Carlo (MC) particle transport methodologies. 

In the post-Apollo era, space operations became somewhat more routine. The increase in 
number of astronauts and the frequency of space visits to Low Earth Orbit (LEO) stations such as 
Skylab, Mir and the International Space Station (ISS) gradually required a reassessment of the 
astronaut protection practices with emphasis on deeper understanding of the potential long term 
genetic mutation and carcinogenic tumor development processes (ref. 4). 

Generally, the RBE for such biological endpoints within the human body, where cancer 
based genetic mutations may initiate, can be large. This is especially true for the types of 
secondary radiation produced in the spacecraft materials and human tissues (refs. 5, 6), pointing 
to the fact that the older dosimetric methods based on tissue equivalent ion chambers (e.g. dose) 
are potentially inadequate. For example, for the absorbed dose in a water shield, figure 1 
compares the analytic fit (nominal 600 MeV protons), derived through the build up factor 
formalism (ref. 6), with experimental (592 MeV proton beam) dose measurements (ref. 7). The 
figure indicates that the formalism of build up which accounts for the nuclear contribution to the 
dose, correlates very well with the measurement of reference 7. In the figure, the noticeable 
sharp rise in the absorbed dose toward the end of beam range is the Bragg peak. This peak occurs 
because the interaction cross section between the beam and the target increases as the beam 
energy decreases through atomic and nuclear collisions. The very sharp rise in the absorbed dose 
at the Bragg peak is often exploited in proton and heavy ion radiotherapy for the treatment of 
cancer. Also shown are the uncollided primaries, where effects of nuclear collisions are ignored. 
From the figure, it is obvious that the secondary particle contribution to the absorbed dose due to 
nuclear processes, especially at low depths, is important in radiation risk estimation for such 
particles and can’t be ignored. 
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depth of water, g/cm 2 

Figure 1 . Absorbed dose in a water shield for a 592 MeV proton beam. 

Ignoring the nuclear processes, the atomic interactions are mainly related to the number 
of electrons per unit volume of material, with small material specific electron binding 
corrections. Consequently, from a practical design point of view, shield properties neglecting 
nuclear reactions can be scaled from one material to another with a fair degree of accuracy. In 
the field of space vehicle design, it was once customary to express all shield materials in terms of 
aluminum equivalent by scaling the effective shield thickness of each material type in terms of 
the range of a 50 MeV proton as follows 

Xeq.Al = [Rai( 50) / R m (50)] . X m (1) 

where Rai(50) is the calculated range of a 50 MeV proton in aluminum (~ 3 g/cm ) and R ra (50) 
in g/cm 2 is the corresponding range in the specific shield material. The ratio in the brackets in 
equation 1 then scales the thickness of material X m in g/cm 2 to an equivalent thickness of 
aluminum X eq .Ai in g/cm 2 , which will have the same shielding effectiveness if nuclear reactions 
are negligible. It is worth noting that by expressing material range in g/cm 2 in equation 1, the 
effect of material density in g/cnT is automatically incorporated in the range calculation. 

The approximation of ignoring nuclear reactions works well for the estimation of 
absorbed dose from impinging protons below roughly 100 MeV, where the proton range is only a 
few g/cm 2 . but begins to generate error for dose at higher energies in the Galactic Cosmic Ray 
(GCR) energy domain. To emphasize the consequence of excluding nuclear reactions in an 
absorbed dose calculation, we return to the 1 60 g/cm water column measurement of reference 7 
as was shown in figure 1 and compute the percent nuclear contribution to the absorbed dose for 
proton beams of 100 ~ 1000 MeV. This is shown in figure 2 where low energy protons have a 
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depth in water, g/cm 2 

Figure 2. Percent dose contribution due to nuclear reactions for a proton beam of varying energy 

impinging on a water column. 

relatively small nuclear contribution to dose and only above 500 MeV or so begin to dominate 
the dose contribution due to nuclear reactions. This implies that for the energetic GCR ions in the 
GeV range, the approximation of equation 1 no longer holds. It is worth noting that the initial 
build up of the absorbed dose is due to secondary neutron production and noticeably for the 
absorbed dose, all nuclear interactions cease to exist near the Bragg peak. 

In the free space environment, where the energies of proton and helium components of 
the GCR reach a few GeV, the inclusion of nuclear processes in any mission planning dose 
estimation becomes more important due to the potential long duration of the mission. Even more 
importantly, in the GCR enviromnent, the problem of estimating radiation health risks to the 
crew of the spacecraft is also strongly particle type dependent, and hence new methods for risk 
estimation must also include the establishment of an accurate particle environment (refs. 8, 9). 
Hence, for any space mission outside the protective geomagnetic field of the earth, great care in 
evaluation of the particle environment is the underpinning of astronaut health risk estimation 
(ref. 10). 

In spite of the advances in computational procedures, due to simplicity of the aluminum 
equivalent approximation, many calculations are still being made using the equivalent 
approximation assumption (ref. 11). This assumption is mainly made out of convenience using 
readily available Monte Carlo (MC) databases generated many years ago. It is, therefore, useful 
to examine what limitations this assumption may have on risk estimate, which is the purpose of 
the present study. 
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2. Computational details 


In this study, the aluminum equivalent approximation is evaluated as a means for 
estimating the particle environment induced by the GCR radiation field within a spacecraft 
structure. For comparison, the two extremes of the GCR environment, the 1977 solar minimum 
and 2001 solar maximum, are considered. Prior work studied these extreme GCR environments 
in detail (ref. 12). 

Models of free space GCR environment developed in the past two decades have provided 
the most realistic description of the interaction of incoming GCR from outside the heliosphere 
with solar activity. The model of O’Neill (ref. 13) is currently used as GCR input into the most 
recent version (2005) of the Langley Research Center (LaRC) deterministic ionized particle 
transport code High charge (Z) and Energy TRaNsport (HZETRN) (ref. 14). This GCR model is 
based on fitting the existing balloon and satellite measured energy spectra from 1954 -1992 and 
more recent measurements from the Advanced Composition Explorer (ACE) satellite from 1997 
- 2002 to the stationary Fokker-Planck equation to solve the diffusion, convection, and energy 
loss boundary value problem and obtain an estimate of the appropriate diffusion coefficient. In 
addition, correlation of the diffusion coefficient to the Climax neutron monitor data which 
exhibits an odd-even cycle with a 22 year period, enables the estimation of the diffusion 
coefficient at times that direct observational data are not available. The latest implementation of 
this GCR model (2004) accurately accounts for the solar modulation of hydrogen through nickel 
(H - Ni) by propagating the Local Interplanetary Spectrum (LIS) of each element through the 
heliosphere. The model provides a single value of the deceleration parameter cp(t) describing the 
level of solar cycle modulation and determines the GCR differential energy spectrum for 
elements 1 < Z < 28 at any given radial distance from the sun. 

The coupling of the GCR extreme environments to HZETRN allows for the examination 
of the induced environment within a spherical shell shield and the effects of the aluminum 
equivalent approximation for a good polymeric shield material such as generic polyethylene 
(PE). For numerical simulation, the shield thickness is restricted to 25 g/cm 2 . Although one 
could imagine the progression to greater thickness, the current range will be sufficient to 
evaluate the qualitative usefulness of the aluminum equivalent approximation. 

Upon establishing the inaccuracies of aluminum equivalent approximation through 
numerical comparisons of the GCR radiation field attenuation through PE and aluminum 
equivalent PE spherical shells, we further present results for a limited set of hydrogen rich 
polymeric constituents such as PE foam core, Spectra fiber composite facesheet, epoxy resin 
block and polyurethane carbon facesheet to assess the effect of the approximation in their 
radiation attenuation response as compared to generic PE. These off the shelf, commercially 
available, multifunctional polymers have promising physical and chemical characteristics as 
potential candidates in space applications. Through experimental testing and numerical 
simulation, prior reports have discussed the structural, thermal and radiation properties of these 
multifunctional polymers in great detail (refs. 15, 16). 
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3. Results 


The transport of the 1977 and 2001 GCR spectra in PE and aluminum equivalent PE 
spherical shell shields was evaluated and percent relative error of particle spectra in PE to that of 
aluminum equivalent PE at the same equivalent thickness for Z=0, sum of 1 < Z < 2, sum of 3 < 
Z < 10, sum of 11 < Z < 20 and sum of 21 < Z < 28 were evaluated. These charge dependent 
relative errors in flux provide an estimate of the limitation of the aluminum equivalent 
approximation for different charge groups in PE. In addition to charge group relative error 
computation, the percent relative error for dose was calculated, as this quantity represents the 
charge and energy integrated limitation of the equivalent approximation. Table 1 represents the 
computational spatial grid for PE and aluminum equivalent PE as described by equation 1. In 
generating table 1, values of R a i(50)=2.93 g/crn 2 and Rpe(50)=2.07 g/cm 2 as calculated by 
HZETRN, were used. This grid was used to propagate ions for both GCR extremes. 


PE (g/cm 2 ) 

Al-eq. PE (g/cm 2 ) 

0 

0 

0.1 

0.14 

0.3 

0.42 

0.7 

0.99 

1 

1.41 

3 

4.22 

5 

7.04 

7 

9.86 

10 

14.08 

15 

21.12 

20 

28.16 

25 

35.2 


Table 1. Computational spatial nodes. 


3.1 Neutron flux error 

The percent relative error for the neutron flux for the 1977 and 2001 GCR solar minimum 
and maximum environments at various shield thicknesses are shown in figure 3. The figure 
indicates that for both extreme enviromnents error below 10 MeV is on the order of 50 percent or 
more with escalating error at lower energies and greater depths. Nearly a factor of two (100 
percent) error occurs at the biologically important energy of 1 MeV beyond 25 g/cm 2 of 
aluminum equivalent PE. At this energy, the proton range in tissue (water) is approximately 
equal to the skin thickness (0.01 g/cm 2 ). Finally, the figure indicates that the flux ratios during 
the 2001 solar maximum are qualitatively similar to those of 1977 solar minimum. 
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Figure 3. Percent relative error for neutron flux for the 1977 and 2001 GCR 
environments (PE vs. aluminum equivalent PE). 


3.2 Light ion flux error 

The corresponding flux percent error for the combined light ions (proton, deuteron, triton, 
helion and helium) is shown in figure 4. In the figure, the error due to the aluminum equivalent 
approximation is on the order of 40 percent at 10 MeV and decline to smaller values at higher 
energies where the penetrating primary ions dominate the spectrum. Below 10 MeV, the error 
rises to over 50 percent at 1 MeV with a continued rise to higher percentages at the lowest 
energies. 



1.0E-02 1.0E-01 1.0E+00 1.0E+01 1.0E+02 1.0E+03 1.0E+04 1.0E+05 

energy, MeV 

Figure 4. Percent relative error for the combined light ion (proton, deuteron, triton, helion and 
helium) flux for the 1977 and 2001 GCR environments (PE vs. aluminum equivalent PE). 
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3.3 Light-medium ion flux error 


Figure 5 shows the combined charges 3 through 10 ion group percent error. The figure 
indicates that the error is on the order of 50 percent at the lowest energies and declines to about 
20 percent at 10 MeV. The figure also indicates that the charge group has a gradual error build 
up at higher depths and energies. The increase in error at higher depths and energies is due to the 
collisional energy losses of the heavier elements in the charge group and is expected to increase 
for even heavier elements. 



1.0E-02 1.0E-01 1.0E+00 1.0E+01 1.0E+02 1.0E+03 1.0E+04 1.0E+05 

energy, MeV 

Figure 5. Percent relative error for the combined charges 3 through 10 ion group flux for the 
1977 and 2001 GCR environments (PE vs. aluminum equivalent PE). 

3.4 Medium and heavy ion flux error 

The relative error for the remaining ion groups, the combined 1 1 through 20 charge group 
and the combined 21 through 28 charge group, are shown in figures 6 and 7, respectively. Their 
qualitative behavior is similar to the charge 3 through 10 group with a noticeable rise in the 
percent error at higher depths and energies. Figure 6 indicates that at a depth of 25 g/cm 2 , the 
percent error across the entire GCR energy range exceeds 30 percent, while in figure 7, at the 
same depth, due to the presence of the even heavier GCR constituents and higher collisional 
losses, the error remains above 50 percent. 

The flux relative error comparison in the previous sections indicates that, with the 
exception of heavier ion groups (11 < Z < 20 and 21 < Z < 28) where collisional losses are large, 
the GCR ions with energies above 1 GeV easily penetrate the maximum depth of the target 
material in this study (25 g/cm 2 ). This was shown by the near zero relative error above 1 GeV for 
lighter charge groups. Flowever, for a large habitable structure (e.g. ISS) where in some 
directions the areal thickness are well above 25 g/cm 2 , this small relative error can become 
considerably larger. 
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energy, MeV 


Figure 6. Percent relative error for the combined charges 1 1 through 20 ion group flux for the 
1977 and 2001 GCR environments (PE vs. aluminum equivalent PE). 
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Figure 7. Percent relative error for the combined charges 21 through 28 ion group flux for the 
1977 and 2001 GCR environments (PE vs. aluminum equivalent PE). 
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3.5 Dose error 


For dose error analysis, the absorbed dose (D) due to the energy deposition at a given 
location depth x of the target material by all GCR ions, is calculated according to 

D(x) = Z j J(p j (x,E)S j (E)dE (2) 

where cpj (x,E) is the flux type j ion at depth x with kinetic energy E, and Sj(E) is the stopping 
power of type j ion with kinetic energy E. 

Applying equations 1 and 2, the percent relative error in absorbed dose using PE and 
aluminum equivalent PE as targets, for various shield thicknesses, is shown in figure 8. The 
figure indicates that using PE and aluminum equivalent PE as targets, the relative error for the 
2001 solar maximum is larger than the 1977 solar minimum. This error is not to be confused 
with the actual absorbed dose values for each epoch where the 1977 solar minimum produces 
larger absolute dose values than the 2001 solar maximum. 



Figure 8. Percent relative error for the absorbed dose for the 1977 and 2001 GCR environments 

(PE vs. aluminum equivalent PE). 

Figure 8 indicates that percent relative error due to the approximation, for both GCR 
extremes integrated across the charge and energy range, is increasing even at the limiting PE 
depth of 25 g/cm 2 of this study. For a large habitable spacecraft module this error can result in 
incorrect absorbed dose estimation at a given target point within the module. 

It is worth noting that the magnitude of the absorbed dose error is significantly less than 
the flux error of previous sections. This is due to the integrated nature of equation 2 where large 
errors at low energies are compensated by small errors at higher energies. 


10 


3.6 Multifunctional polymer error 


In previous studies (refs. 15, 16), HZETRN, in conjunction with appropriate 
environments (ref. 12), was used to explore the response of commercially available, hydrogen 
rich polymers to ionizing radiation, as a first step toward the long term goal of developing 
multifunctional composite structural configurations and manufacturing methodologies to reduce 
spacecraft weight while maintaining structural, damage tolerance and radiation shielding 
requirements. The best such hydrogen rich multifunctional polymers with acceptable structural 
properties are aliphatic polymers such as high-density PE fibers, epoxy matrix/adhesive and PE 
foams. Such polymers will most likely have to be manufactured in the form of sandwich 
structures where the lightweight and radiation tolerant core is supported by stiff fiberglass/epoxy 
or carbon/epoxy facesheets to form a structurally efficient membrane. 

In this section, we apply the approximation to some of these commercial grade polymers 
to develop a quantitative understanding of the extent of relative error in the absorbed dose. We 
will then compare these errors with the PE error of the previous section to draw conclusions 
about whether the approximation should be used for polymers in general. 

Table 2 is an abbreviated list of the important physical and chemical properties of the 
seven candidate polymers studied in reference 16. In the table, elemental percentages in red refer 
to the highest percentile concentration of elements H. C, N and O. 


Material 

Description 

Density 

(g/cm 3 ) 

Chemical 

Composition (mass%) 

Polyethylene 

generic PE 

0.94 


Resin 1 

Epoxy resin block 

1.12 

70.1% C 7.6% H 16.8% O 5.5% N 

Facesheet 1 

Spectra/epoxy facesheet 

1.04 

77.3% C 10.7% H 9.0% O 3.0% N 

Facesheet 2 

Carbon/epoxy composite 

1.44 

83.5% C 3.6% H 7.10% O 5.8% N 

Core 1 

PET* foam block 

0.35 

62.5% C 4.2% H 33.3% O 

Core 2 

PP* foam block 

0.19 

85.6% C 14.4% H 

Sandwich 
Composite 1 

0.35g/cm 3 PET*/Spectra 
facesheet 

0.52 

69.9% C 7.4% H 21.4% O 1.3% N 

Sandwich 
Composite 2 

FR6720 Polyurethane/carbon 
facesheet 

0.53 

75.2% C 4.6% H 13.2% O 6.4% N 
0.3% P <1% Trace elements 


PET : Polyethylene Terephthalate 

* 

PP : Polypropylene 


Table 2. Properties of the target multifunctional polymers investigated. 
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Since the aluminum equivalent approximation is based on the 50 MeV proton range in 
any material, table 3 provides the 50 MeV proton range in all the materials listed in table 2. In 
addition, the range in aluminum is also included as a reference. The ranges in the table were 
calculated using HZETRN. For PE and aluminum, the ranges were checked versus National 
Institute of Standards and Technology (NIST) predictions (ref. 17). 


Material 

50 MeV proton range 
(g/cm 2 ) 

Polyethylene 

2.07 

Aluminum 

2.93 

Resin 1 

2.25 

Facesheet 1 

2.16 

Facesheet 2 

2.37 

Core 1 

2.37 

Core 2 

2.07 

Sandwich Composite 1 

2.26 

Sandwich Composite 2 

2.35 


Table 3. 50 MeV proton range for the target multifunctional polymers investigated. 

With properties of these multifunctional polymers defined, we now calculate the percent 
relative error in absorbed dose for the 1977 and 2001 GCR spectra. Figure 9 displays percent 
relative error for the absorbed dose. The largest errors for both epochs are due to PE and Core 2 
polymers which have very similar properties. Note that the percent absorbed dose relative error 
of the 2001 solar maximum is somewhat larger than the 1977 solar minimum. As discussed 
previously, this error is not to be confused with the actual absorbed dose values for each epoch 
where the 1977 solar minimum produces larger absolute dose values than the 2001 solar 
maximum. 

Figure 9, which in essence is a composite representation of aluminum equivalent 
approximation to most polymers, indicates that the application of this approximation to polymers 
results in incorrect dose calculation and should be avoided. Note that in the figure, even at a 
depth of 25 g/cm , which is most likely well below the 50 percentile average thickness of ISS for 
the 2008 configuration of 1J (ref. 18), the error is about 16 percent and increasing. This error 
level points to the problem that use of the equivalent approximation can result in incorrect dose 
estimation at a given target point. 
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1977 GCR solar minimum 


2001 GCR solar maximum 
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— 

PE 

2 

— 

Resin 1 

3 

— 
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— 
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5 

— 
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Sandwich Composite 1 

8 

— 

Sandwich Composite 2 


Figure 9. Percent relative error for the absorbed dose for the 1977 and 2001 GCR environments 
(multifunctional polymer vs. aluminum equivalent multifunction polymer). 


4. Conclusion 

This paper elaborated on the historical justifications for the sole emphasis on atomic 
interactions, and why nuclear collisions were essentially ignored. Also discussed was the fact 
that from a vehicle design point of view, shield properties, neglecting nuclear reactions, can be 
scaled from one material to another with a fair degree of accuracy, relative ease, and huge saving 
on computational time. The relative ease of material scaling combined with the acceptable turn 
around time for numerical simulation were indeed the impetus behind the equivalent 
approximation concept for the first thirty years of the space program. 

The aluminum equivalent shield approximation was a matter of convenience when dose 
was the monitored quantity to guard against potential radiation sickness in the early space 
program. With soon to be established permanent human presence in space (ISS) and potentially 
a lunar base, and the search for mitigation methodologies to minimize exposure by appropriate 
choice of shielding materials, applying the equivalent approximation to polymers is no longer 
justifiable. Under the GCR exposures considered herein, this study showed that for polymer 
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based shielding materials, the approximation produced large errors in flux and absorbed dose 
estimations. 


Clearly, with the advent of fast numerical marching techniques such as HZETRN and 
other deterministic algorithms and the rapid maturation of high performance computers, the use 
of this approximation is unjustifiable for polymers. 
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